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ABSTRACT 

The ionized core in the Sgr B2 Main star-forming region was imaged using the Submillimeter Array 
archival data observed for the H26a line and continuum emission at 0.86 millimeter with an angular 
resolution 0.3". Eight hyper-compact H26a emission sources were detected with a typical size in 
the range of 1.6-20xl0 2 AU and electron density of 0.3-3xl0 7 cm -3 , corresponding to the emission 
measure 0.4-8.4xl0 10 cm -6 pc. The H26a line fluxes from the eight hyper-compact HII sources 
imply that the ionization for each of the sources must be powered by a Lyman continuum flux from 
an O star or a cluster of B stars. The most luminous H26a source among the eight detected requires 
an 06 star that appears to be embedded in the ultra-compact HII region F3. In addition, ~ 23 
compact continuum emission sources were also detected within the central 5"x3" (~ 0.2 pc) region. 
In the assumption of a power-law distribution for the dust temperature, with the observed brightness 
temperature of the dust emission we determined the physical properties of the submillimeter emission 
sources showing that the molecular densities are in the range of 1-10 xlO 8 cm -3 , surface densities 
between 13 to 150 g cm" 2 , and total gas masses in the range from 5 to > 200 Mq which are 1 or 2 
orders of magnitude greater than the corresponding values of the Bonnor-Ebert mass. With a mean 
free-fall time scale of 2xl0 3 y, each of the massive protostellar cores are undergoing gravitational 
collapse to form new massive stars in the Sgr B2 Main core. 

Subject headings: Galaxy: center — ISM: individual (Sgr B2) — HII region: radio lines — radio 
continuum: ISM — stars: formation 



1. INTRODUCTION 

The formation of massive stars is one of the chal- 
lenging problems in astrophysics. Unlike their low- 
mass countparts, massive stars are rare, and form in 
relatively deeply embedded m assive molecular clouds 
on a much shorter time-scale { McLaug hlin k Pu dritz 
19971 lOsorio. Lizano fc D'Alessio 11199a iMcKee fe Tan I 

20021 ) . High-mass sta rs are often formed in 

clusters (|McKee k Ostrikerl 120071 ). After initial grav- 
itational collapse of the natal clouds, multiple proto- 
stellar cores are formed through fragmentation of the 
gas. One of the fundamental theoretical questions is 
how the massive stars f orm in clusters. Accordin 
to the competitive theory (IBonnell. Vine, fc Bate 11200 
iBonnell fc Bate 1120061 iBonnell 1120081 ). mass clumps cre- 
ated from natal cloud gravitational collapse contain pro- 
tostellar cores with small masses. These cores grow 
by accreting matter, competing with other cores, and 
form stars with many times their original mass. On 
the other hand, the direct gravitati onal collapse the- 
ory (jKrumholz. McKee k K lein 2005) suggests that the 
protostellar cores created from fragmentation of the ini- 
tial cloud gravitational collapse have sufficient mass to 
form individual high-mass to low-mass star systems in 
the subsequent collapse. Accretion from the parental 
cloud continues but does not substantially change their 
mass. 

Strong radiation pressure from a newly formed mas- 
sive star might halt infall, limiting the mass of stars 
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that can form (IKahn lH9jl IWolfire fc Cassinelli I U987I: 
ILarson fc Starrfield 1119711 7 Recent theoretical investiga- 
tions suggest that the radiation pressure limit might be 
overcome in actual cases with complex, non-spherical in- 
fall geometries or hig h ram pressures in rota t ing disks 
(IMcKee k Tan 1 120021: iKrumholz et al. I l200a iNakano I 
1989), or massive st ars may form from stellar merg- 
ing in a dense cluster (IBonnell. Bate, k Zinnecker 1119981 : 
22). 



IBonnell fc Bate 2002). Numerical simulations suggest 
that graviational instabilities cause the disk to frag- 
ment and form a mas sive companion to the primary 
( Krumhol z et al. Il2009ft , a nd consequently form b inaries 
in a dense stellar system (|Bonnell fc Bate 1 12002ft . Fur- 
thermore, radiation feedback from massive stars in a clus- 
ter affects its subsequent fragmentation and consequently 
plays an important role in determini ng the stellar ini- 
tial m ass function (IMF) in a cluster (jKrumholz et al.l 
[20101) . 

In the past decade, good progress has been made in un- 
derstanding mas sive young st ellar clusters in the Galaxy. 
A recent review (|Figer 1120081 ) shows ten known Galactic 
clusters with masses < 10 4 Mq with ages of a few million 
years. The Arches cluster in the Galactic center is the 
densest young cluster in t he Galaxy an d contains a large 
collection of massive stars (jFiger II20"05T ) . Three out of ten 
(Quintuplet, Arches and Center) are located within the 
central 50 pc of the Galaxy, suggesting that the Galactic 
central region appears to prefer forming ma ssive stars. 

Lo cated at a distance of 7.8 kpc ()Reid et al. I 
2009) towards the Galactic center, Sgr B2, a gi- 
ant molecular cloud(GM C) with a mass of 6xl0 6 M Q 
(|Goldsmith fc Lis I [l990h is one of the most active 
star-forming region in the Galaxy, radiating a to- 
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tal luminosity o f lxl 7 L w (|Goldsmith fc Lis I 119901: 
iGoldsmith et al. I Q.992). As the most luminous core 
among the several cores in this GMC, Sgr B2 Main 
is associated with numerous ultra-compact (UC) HII 
regions, s uggesting the presence of a tight cluster of 
OB stars (IGaume fc ^ Claussenlll990l IGaume et "all 119951; 
iDe Pree et al. Ill998f) . The observed molecular outflows 
and infalls sugge st that on-going, star formation acti vities 
are taking place (|Qin et al.ll2008l: iRolffs. et al.ll2010[) . Sgr 
B2 Main appears to be in a very young phase of forming 
a massive stellar cluster from the dense molecular core. 

High-resolution observations using the Submillimeter 
Array (SMA j3 at submillimeter wavelengths can explore 
the detailed structure of the massive star forming core, 
providing useful clues on how massive stars form in a 
cluster. 

The rest of this paper is organized as follows. In Sec- 
tion 2, we discuss the reduction and imaging process of 
the SMA archival data of Sgr B2 Main observed in 2007. 
Section 3 shows the results from the high-resolution ob- 
servations of the H26a line and continuum at 0.86 mm. 
Section 4 presents a model for the H26a sources. In Sec- 
tion 5, the properties of the bright, compact dust emis- 
sion clumps are determined and discussed. In Section 
6, we discuss the kinematics in terms of ionized out- 
flows/expansions and rotating disks. The early phase 
of massive star formation, and origin of the massive pro- 
tostellar clumps are also discussed. In Section 7, we sum- 
marize the conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

The interferometer data for the H26a line at vq — 
353.623 GHz were acquired from the SMA archive, 
observed 2007 June 18, with the "very extended" 
array configuration in the upper-side band (USB). 
The reduction for the line dat a was made in Miriad 
(jSault. Teuben. fc Wrightll995l) following the reduction 
instructions for SMA dat£0. The bandpass calibra- 
tion was made by applying a normalized average of 
all the data from the QSOs included in the observing 
run (J1229+020, J1733-130, J1743-038, J1751+096 and 
J2015+371). The flux density scale was determined by 
comparing observations of Callisto with the SMA plane- 
tary model. An angular size and brightness temperature 
of 1.5" and 120 K at the observing epoch were assumed. 
The complex gains of the data were primarily calibrated 
using the nearby QSO J1733-130 to obtain an image of 
the SgrB2 Main core. This image was used as an initial 
model to further correct residual errors in the visibilities 
using self-calibration. 

To separate line and continuum emission, we used the 
Miriad task UVLIN. The line visibility data set was con- 
structed by subtructing the continuum which was deter- 
mined from line-free channels. We made images of the 
H26a line with a channel width 3 km s -1 using a robust- 
ness weighting of 2 corresponding to natural weighting. 
The FWHM beam is 0.40" x 0.28" (PA=15°). The typ- 
ical rms noise in a channel image is 0.1 Jy beam -1 . 

3 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 

4 http://www. cfa. harvard, edu/sma/miriad 



A line-free continuum data set was also produced from 
the UVLIN program. We used both USB and LSB data 
to image the continuum emission with robustness weight- 
ing parameters 2 and -2 which correspond to natural and 
uniform weighting respectively. The FWHM beam with 
uniform weighting is 0.36" x 0.22" (PA=13°). The typ- 
ical rms noise in the continuum images are 7.5 and 8.5 
mjy beam -1 for natural and uniform weighting, respec- 
tively. The shortest sampled spatial frequency in the 
visibility data is 30 kA, corresponding to an angular size 
~ 6"-7" for the H26a line and continuum emission struc- 
ture. The compact structure (<3") in the H26a line and 
continuum emissions from the Sgr B2 Main core has been 
adequately sampled in the SMA observations. 

We noticed a systematic offset of 0.3" in the phase 
center of the SMA images by comparing the positions 
of the UC HII regi ons G and F10.37 d etermined at 22.4 
GHz with the VL A (jGaume et al.ll995| ) . The offset could 
reflect residual errors in antenna positions and complex 
gain calibration due to the large distance between the 
calibration QSO J1733-130 and SgrB2 Main. We there- 
fore imposed a shift of 0.3" for the phase center of the 
SMA data to align with the VLA image. The positional 
error of a compact source with S/N — 6 in the final 
SMA images with respect to the VLA coordinate frame 
is -0.03" (see below). 

3. RESULTS 

3.1. H26a line emission 

3.1.1. Distribution 

Figure QJi shows the integrated H26a line image of Sgr 
B2 Main made from line emission in the LSR velocity 
range between —100 and +225 km s -1 . Eight sources 
with significant H26a line emission (> 6cr in the chan- 
nel image) are labelled as H26a-n, where n is numbered 
from 1 to 8. The source H26a-1 has a small angular offset 
from the brightest continuum source (F3) at 22.4 GHz 
in the core of Sgr B2 Ma in while H26a-8 ha s a larger 
angular distance from F3 (jGaume et al.l fl995). The two 
relatively isolated bright line sources, H26a-8 and H26a- 
5, were used to align the 0.86 mm coordinate frame of our 
SMA images with the coordinate frame dertermined with 
the VLA at 22.4 GHz. The image shows the submillime- 
ter positions of H26a-8 and H26a-5 agree with the 22.4 
GHz positions of G and F10.37 within -0.01". H26a- 
6 appears to coincide with F10.303. However, H26a-1, 
the brightest line source, is located 0.17" northwest of 
F3. H26a-6 coincides with F10.303. Two ultra-compact 
H26a line sources (3 and 4) appear to be associated 
with Fl. The hyper-compact H26a sources 3 and 4 
are located —0.15" northeast and southwest of Fl, re- 
spectively. Also, H26a-2 is located 0.2" southwest of 
F4. No significant H26a emission associated with F2 has 
been detected. The positional offsets between the inner 
four H26a emission sources (1,2,3 and 4) and their corre- 
sponding 22.4 GHz counterparters (F3, F4 and Fl) ap- 
pear to be signficant, reflecting the fact that the free- free 
continuum emission at 22.4 GHz traces relatively lower- 
density ionized gas in the outer layer of the expanding 
HII gas or the ionized outlows, while the H26a line traces 
the high-density ionized regions, possibly associated with 
the ionized disks or ultra- or hyper-compact ionized cores 
where the free-free emission appears to be optically thick 



SAGITTARIUS B2 MAIN 



3 



H26a-8 




L - 







H26a-2 




H26a-5 




H26a-7 





H26a-6 



H26a-3 f 







1 00 




-109 



28 23 08 



17 h 47 m 20 s 3 


20 s 2 20 s l 


20 s 




RA (J2000) 




r 1 i i i i | i i i 
I i i 



H26a-1 I 







30 



Jy/beam km/s 



100 



150 



Fig. 1. — a: SMA image of H26« line emission at the rest frequency of 353.623 GHz, integrated from the velocity range of —100 to 
125 km s — 1 (red contours and color) overlaid on the continuum image at 0.86 mm (yellow contours). The red contours are i<p(H26a) X (4, 
5, 7, 11, 17, 25, 35, 47, 61, 77 and 95%), and the yellow contours are S P (0.86mm)x (-2, 2, 3, 5, 10, 15, 25, 35, 45, 55, 65, 75, 85, and 
95%), where the peak line flux i 7 p(H26a) and peak continuum flux density Sp (0.86mm) are 193 Jybeam -1 kms — 1 and 2.88 Jybeam -1 , 
respectively. The FWHM beam is 0.40" xO. 28" (15°) denoted at bottom-right. The color wedge scales the H26a line intensity in the units 
of Jybeam - kms" 1 of the UC HII regions. The spectra from the eight brightest H26« sources are also shown. The units of the vertical 
(line flux density) and horizontal (LSR velocity) axes in the spectral plots are Jy and kms -1 , respectively. The red cross in each spectral 
profile shows the peak line intensity (vertical bar) and the FWHM velocity width of a Gaussian profile determined from least-square fitting. 

at 22.4 GHz. iBrocklehurst k Seatonl Il972f ). We find that the pres- 

sure broadening (AVp ~ 0.02 km s -1 ) is negligible for 
the H26a line. Thus, the non-thermal broadening AV nt h 
including the Doppler motions of expansion, infall, out- 
flows, rotation, shocks, and turbulence, obtained by sub- 
tracting the thermal broadening from the observed H26a 
line width, is less than 34 km s _1 . 

For the H66a line, the pressure broadening for gas with 
n c = 1.1 x 10 6 cm -3 accounts for 23 km s _1 . Subtract- 
ing AT4h, AV„th and AVp from the observed AVfwhm 
for the H66a line, ^AV F ' wliM - (AVg + AV^ th + AV^) 
ps 43 km s 1 , the residual line width, appears still to 
dominate the observed H66a line width, suggesting that 
the el ectron density estima ted from the 7 mm continuum 
data ()De Pree et al. |[l~998f) might be underestimated due 
to the large optical depth of the free-free emission from 
the higher density gas. If the residual line width is all 
due to pressure broadening, the observed large H66a line 
width implies that a hyper-compact ionized core with an 
electron density rt c >2.5 x 10 6 cm -3 is likely embedded 
in F3. 



3.1.2. Line broadening 

The spectral profiles of the H26a line from the eight 
brightest UC HII regions in the core of Sgr B2 Main 
are shown in Figure [TJ H26a-1, associated with F3, 
has a peak velocity Vlsr = 64 ± 1 km s -1 with a 
line width AVfwhm = 40 ± 1 km s _1 compared with 
Vlsr = 68 ± 2 km s _1 and AVfwhm = 63 ± 5 km 
s _1 derived from the VLA measurements of the H66a 
line (|De Pree et al. 1119961) . The peak velocity between 
the H26a and H66a line profiles shows a 2a difference. 
The line width of the H66a is much broader than that 
of H26a. Using canonical values of T c = 1 x 10 4 K 
and n c = 1.1 x 10 6 cm -3 for the typical U C HII re- 
gions in Sgr B2 Main (|De Pree et al 1 119981 ). for the 
VLA observations of the continuum emission at 7 mm, 
the thermal broadening AVth contributes 21 km s _1 
to the line widths. The pressure broadening depends 
upon the principal quant um number N and electron den- 
sity n e (AVp oc n e N 7A ) (|Brocklehurst fc Leeman H1971I : 
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Fig. 2. — The image of 0.86-mm continuum emission observed with the SMA at a resolution of 0.36" X 0.22" (PA=13°). The contours 
are (-0.02,0.02,0.05,0.15,0.25,0.35,0.45,0.55,0.65,0.75,0.85,0. 95)x2.34 Jy beam" 1 . The color wedge shows the intensity range in units of Jy 
beam -1 . The positions (cross) and the apparent shapes (ellipses) of the continuum sources are derived from Gaussian fitting. 

peak brightness and total flux densites (columns 5 and 
6). The brightness temperatures of the sources are given 
in column 8. We note that the two brightest continuum 
cores are Smm-1 and Smm-4 with brightness temper- 
atures 340 and 270 K, respectively. Smm-1 is located 
close to but with a siganificant offset 0.13"±0.01" from 
H26a-1 which is associated with F3. Smm-4 is located 
0.06"±0.02" southeast of its H26a counterpart H26a-4 
in the Fl complex. 

4. MODEL FOR UC HII REGIONS 

On the assumption of isothermal homogeneous HII gas, 
a model has been developed for the low-frequency radio 
recombination lines (RRL) including non-LTE effects, 
such as stimul ated emission by the background radiation 
(jShaver II1975T ) . This model has been applied successfully 
to the high-frequency radio recombination lines (H30a) 
at 231.9 GHz from the high-density gas in the minispi- 
ral within the cir cum-nuclear disk of the Galactic center 
(|Zhao et al.ll2010f ) . The effects of the radiation from both 
synchrontron and dust sources seem to be negligible for 



3.2. Continuum Emission 

The contours in Figured] show the continuum emission 
with the same FWHM beam (0.4" x 0.28") as the H26a 
line. The continuum flux densities corresponding to the 
H26a line emission regions are summarized in column 8 
of Table 2. A large fraction of the continuum emission at 
0.86 mm appears to arise from a region outside the H26a 
line sources, revealing an overall extended filamentary 
structure northwest to southeast. Figure [2] shows the de- 
tails of the continuum emission at the higher resolution 
(0.36" x 0.22"). In addition to the UC HII regions ob- 
served at millimeter and centimeter wavelengths at the 
positions denoted with the red dots, we identified 23 con- 
tinuum sources at submillimeter wavelengths which are 
listed in the bottom section of Table 2 and are named 
as Smm- with a sequential number (column 1) based on 
their angular distance Ad from F3 (column 7). The cor- 
responding sequential numbers of the Smm- sources are 
also marked in Figure [2j Using Gaussian fitting, we de- 
termined their positional offsets from F3 in RA and DEC 
(columns 2 and 3), deconvolved sizes (column 4) and 
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the H30a line under the conditions within the central few 
parsecs of the Galaxy. In the massive star forming cores 
such as Sgr B2 Main, a considerable fraction of the con- 
tinuum flux density at submillimeter wavelengths is from 
dust radiation that needs to be included in the model. 
Using a homogeneous isothermal model with tempera- 
tures of T e and Td for both HII gas and dust, the flux 
densities of hydrogen recombination line (Sl) and contin- 
uum (Sq) emission at radio to submillimeter wavelengths 
can be formulated assuming that the total dust is evenly 
distributed in three regions (each with an optical depth 
Td) along the line-of-sight — one third of the dust uni- 
formly mixed with the ionized gas in a given HII region, 
and other two thirds are located in front of and behind 
the HII region, respectively. Equations (1) and (2) give 
the full solutions for the line and continuum flux densi- 
ties from the three regions with no assumptions of small 
optical depths. Spectral fits to the data are shown in 
Figure [3]. The results are not significantly different from 
an optically-thin dust approximation. This is because, 
in the centimeter to submillimeter wavelength range dis- 
cussed in this paper, the dust emission towards the UC 
HII regions is indeed optically thin and we are not able 
to discern between the two models. However, at mid- 
IR (Spitzer) wavelengths, the dust becomes opaque and 
the dust attenuation towards the HII regions becomes 
important. 



Sl — — T c 



tl//?n + t c 
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e~ T c 



B(T d )(l-e- Td )fi 



-7t_! 



<T C +T d )^ (1) 
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1 - e~ T c 
B{T d ) Td n -, e 



~c 

Td 



+ 



B(T d )(l - e~ T *)n 



l + e -(Tc+T d ) 



(2) 



where k is Boltzmann's constant, c is the speed of light, 
and the total continuum optical depth in the HII re- 
gion t'q = tq + Td includes both the free-free and dust 
contributions. The recombination line (tl) and free-free 
continuum (tc) optical depths of the radiation from an 
HII region are given by Equations (A3) and (A4) in Ap- 
pendix A, respectively. The values of the population 
departure coefficients (& n and /3n) were calcula t ed us - 
ing the non-LTE code of iGordon fc Sorochenkol (|2002l ) 
based on the analysis of hydrogen recombination lines 



at wa velengths from radio to submillimetcrs (Walmsley 
1 19901) . The quantity B(T d ) is the Planck function with 
a dust temperature T d . In Equation (1), the first term 

associated with V f flT c is the line emission from the 

c 

HII gas attnuated by the foreground dust with optical 
depth tj; the second term associated with B(T d )T d £l ac- 
counts for the absorption of the internal dust radiation 
by the HII gas if tl > 0, or the line emssion stimulated 
by the internal dust radiation if tl < 0; and the third 
term associated with B(T d )(l — e~ Td )£l is the absorption 
of background dust radiation by the HII gas if tl > 0, 
or the line emssion stimulated by the background dust 
radiation if tl < 0. On the other hand, for an isolated 
dust source, the dust continuum flux density is 



S d = B(T d )(l 



-T d 



)fi. 



(3) 



The optical depth Td of dust continuum radiation can be 
described as 



T d = Ko {-) M d D- 2 n-\ 

v <3 



(4) 



where n v = kq{v/vo)^ is the dust opacity per unit dust 
mass, M d is dust mass, D is the distance, and B v (T d ) is 
the Planck function with a dust temperature of Td- In the 
calculations throughout this paper, we adopted ko = 1-06 
cm 2 g _1 at the reference frequen c y vp = 230 GHz calcu- 
lated by lOsscnkop f fe Henning I (|1994l ) for high-density 
gas riH = 10 7 cm" 3 using the st andard MRN model 
(jMathis. Rumpl fe Nordsiek I I"l977[ ) with thin ice man- 
tles. The physical parameters for a given HII region can 
be determined by fitting an isothermal, homogeneous 
model to the data at the submillimeter wavelengths of 
this paper and the data taken from previous observa- 
tions at radio and millimeter wavelengths as discussed in 
the following section. 

4.1. Radio continuum spectrum & RRL 

F3 is the brightest UC HII region in the Sgr B2 
Main co re as observed at wavelengths of 2, 1.3 and 
0.7 cm (IGaume fc Claussenl 119901 : IGaume et all 119951 : 
iDe Pree et al. II1998D . F3. which is associated with the 
brightest H26a line source, H26a-1, was detected at 2 cm 
with a peak flux density 0.138 Jy beam -1 (with a circular 
beam of 0.3" ), and a total flux dens ity 0.48 Jy in an ex- 
tende d region of size 0.7", FWHM ([G aume & Claussen 
fl990l K With a size 0.5"x0.5", FWHM, the peak in- 
tensity and total fl ux density at 1.3 cm are 0.22 Jy 
beam" 1 and 0.9 Jy (|Gaume et al.lll995l ). Five HII com- 
ponenents have been detected at high-an gular resolution 
(0.04 9"x0.079") with the VLA at 7 mm (|De Pree et ai~l 
1998). The authors show that the major component 
F3-d has a shell- or ring-like morphology. Observa- 
tions with the Hat-Cr eek interferometer array at 3.5 mm, 
iLiu fc Snyder I (|1999f ) showed that the peak emission of 
1.51 Jy beam" 1 (FWHM = l"x0.5") is coincident with 
F3 corresponding to a total flux density of 4.38 Jy for Fl, 
F2, F3 and F4. The high-resolution (0.64" x 0.31") SMA 
image at 1.3 mm shows that the peak emission of 1.8 Jy 
coincides with F3. The H66a line has been detected 
with l.T'xl.6" resolution at V LS R = 68 ± 2 km s -1 
showing a line-to-continuum (L/C) ratio of 3±0.3 % and 
AVfwhm = 63 ± 5 km s _1 determined from integrated 
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Fig. 3.— Models for the UC HII regions F3, F4, F1A, FIB, F10.37 F10.303, F10.39 and G with the constraints from the observed radio 
continuum flux densities and the fluxes of hydrogen recombination lines at radio and submillimeter wavelengths. The red circles denote the 
integrated flux (1.064 S p * AVfwhm) °f the hydrogen reco mbination lines at H26a (from the SMA measurements of this paper), H52q and 
H66« from the VLA measurements (Dc Prcc ct al. 1996). The black dots are the flux densities determined at 2, 1.3 cm and 7 mm with 
the VLA and 3 mm with the Hat Creek interferometer array and 1.3 and 0.86 mm with the SMA (see Section 4.1, the 1.3 mm data from 
Zhu 2010, personal communication). The red curves indicate the fluxes of the hydrogen recombination lines from the best fitted model (see 
section 4) to the observed data. The black curves indicate the continuum spectral energy distribution (SED) resultant from the free-free 
(flat dashed curves corresponding to the first term in Equation (2) ) and dust (steep rising curves corresponding to the second term in 
Equation (2)) emissions calculated from the best fitted model. 



line profile (|De Pree et al. I 119961 ) . At the lower reso- 
lution 2.9"xl.6", the H52a line is towards the F clus- 
ter at Vlsr = 57 ± 2 km s" 1 with L/C of 14±2 % 
and AVfwhm = 59 ± 2 km s^ ljDe Pree et al. I fl996h . 
iDe Pree et al~l (|1996Q noticed that the LTE electron 
temperature (T*) derived from the high L/C value of 
the H52a line is significantly lower by a factor of 2 than 
that derived from the H66a line under the assumption of 
optically thin and LTE gas. If the L/C ratios are uniform 
across the F complex, the fluxes of the H66a and H52a 
lines from the dominant source F3 can be estimated from 
the L/C values and the continuum flux densities at the 
corresponding wavelengths. 

For the other seven regions associated with the H26a 
line source, the continuum flux de nsities at 2, 1.3 and 
0.7 cm used in this paper are fr om iGaume fc Claussenl 
(pQTlh . IGaume et alJ (fl995l ) and IDe Pree et al. I (j!99iC 
The H66a line flux from the UC HII region G (H26a- 
8) is estimate d from the L/C ratio determined by 
IDe Pree et al. I (|1996f) and the continuum flux density 



of IGaume etaLl (fl995h . 

Figure [3] shows the observed spectra of the continuum 
(black dots) and hydrogen recombination lines (red dots) 
at radio and submillimeter wavelengths. 

In order to fit to both the continuum and hydrogen re- 
combination line data obtained from the high-resolution 
observations at wavelengths in the range between 2 cm 
to 0.86 mm, we considered a few possible models. A 
model with an LTE approach and a single isothermal, 
homogeneous HII component was rejected since the ob- 
served large H26a line flux requires large free-free con- 
tinuum optical depths (re > 1) at longer millimeter and 
centimeter wavelengths for a single component model. 
At these wavelengths, the dust absorption (rd << 1) is 
negligible. Then, the predicted line fluxes at centimeter 
wavelengths for larger quantum-number transitions are 
much smaller than the observed values because of the 
exponential drop in the LTE line flux due to the large 

2kv 2 
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The observed line flux densities at centimeter wave- 
lengths could be fitted with a single isothermal, homoge- 
neous component by adjusting n c and T c if the stimulat- 
ing effect in a non-LTE approach is considered. How- 
ever, the fiat continuum spectra at centimeter wave- 
lengths seen in all the eight HII sources (Figure [3]) place 
a critical restriction on a model with a single isother- 
mal, homogeneous component, requiring an additional 
lower-density component. Thus, we fit each of the eight 
observed hyper-compact H26a sources and the UC HII 
regions surrounding them with two components corre- 
sponding to high- and low-density ionized gas with small 
and large sizes, A and B, respectively. Further assump- 
tions used in the two-component model are the volume 
filling factors fy = 0.1 and 0.5 for components, A and 
B, respectively. The effect of shadowing between the two 
components is negligible. 

4.2. Hyper- compact HII component 

Figure [3] shows the results of model fitting the observed 
hydrogen recombination line (red curves) and contin- 
uum flux densities (black curves) from each of the eight 
UC HII regions with two isothermal, homogeneous HII 
components A and B. The best fitted parameters are 
summarized in Table 3. The first row gives the size in 
units of 10 3 AU. The following rows under the category 
of HII gas properties summarize the derived parame- 
ters including the temperature T e , density n c , volume 
filling factor fy, the departure coefficients &n and /3n, 
the optical depths of the H26a line T^(H26a), the free- 
free continuum tq (0.86 mm) at 0.86 mm, the fractional 
contribution Ch26q to the observed H26a line flux from 
each component, the H66a line tl(H66q;), and the free- 
free continuum rc(A13mm) at 13 mm. Components A 
have small size (160 to 2000 AU), high electron density 
(3 to 33 xlO 6 cm -3 ), relatively lower temperature (5 
to 9xl0 3 K). The A components are typically 10 times 
smaller but 100 times de nser than typical U C HII regions, 
10 17 cm and 10 4 cm -3 , (|Churchw ell 2002j), representing 
a cla ss of hyper-compact HII components (|Churchwelll 
2002) present in the Sgr B2 Main core. The hyper- 
compact HII components could arise from plausible ion- 
ized disks with emission measure (EM) in the range 0.4 
to 8.4xl0 10 cm -6 pc, accounting for the most of the 
H26a line flux, e.g. ^94% for F3. The hyper-compact 
HII components in all the HII cores show a negative op- 
tical depth for the H26a line, suggesting the presence of 
stimulated emission by free- free continuum emission from 
the hyper-compact HII region itself. For F3, the H26a 
line is enhanced by a factor of 1.4 due to a weak stimu- 
lating process within the hyper-compact HII component 
(A) while the enhancement factors (Sl/Sl) are 1-^ ana - ^ 
for the H26a lines from G and F10.39 which show larger 
negative line optical depths of —0.71 and —1.3, respec- 
tively. At 13 mm, the continuum component A becomes 
optically thick while the H66a line is still optically thin, 
T C (13mm) = 9.5 and r L (H66a) = 0.13 for F3. Compo- 
nent A appears to make a little contribution to the total 
line flux of the H66a transition (~20% for F3). 

4.3. Ultra-compact HII component 

The continuum emission at longer wavelengths from 
each of the UC HII regions shows a shallow rising spec- 



tra towards short wavelengths, indicating the presence 
of a larger HII component which might result from an 
ionized stellar wind or expansion of an ionized nebula. 
Component (B) with a relatively large size (1200 to 4000 
AU), lower density (0.1 to 3.5xl0 5 cm -3 ), and higher 
temperature (10 to 17x 10 3 K) is needed for the extended 
ionized gas with relatively smaller emission measures of 
< 1 x 10 9 cm- 6 pc. T c = 17, 000 K for the UC HII 
component in F3 (H26a-1) appears too high for the HII 
regions with metallicities in the Galactic center. The 
over-estimated T c from our fitting might be caused by 
the under- estimate of the H66a and H52a line fluxes at 
1.3 and 0.7 cm due to large line broadening and lim- 
ited bandwidth coverage of the old VLA system. For 
a canonical value of T c — 10, 000 K, the corresponding 
H66a and H52a line fluxes require 50% more than the 
values used in this paper. The physical parameters of the 
B components fal l into the category of UC HII region in 
iChurchwelll (|2002ft . The continuum emission from the B 
components is optically thin at longer millimeter wave- 
lengths, providing a flat spectrum (oc v~ ) in contrast 
to the steep rising spectrum (oc u ) of component A. 
With its larger size, the free-free emission from the lower 
density gas (component B) dominates the continuum flux 
densities at 13 mm. Superposition of the spectra from 
the two components (A and B) results in a spectrum 
slowly rising at long wavelengths and turning over to flat 
at short millimeter and sub-millimter wavelengths. 

4.4. Dust component 

The composite spectra (dashed flat curves in Figure [3]) 
of components A and B appear to fit the observed con- 
tinuum flux densities at millimeter and centimeter wave- 
lengths. In the submillimeter, the observed flux densities 
show a significant excess emission from the UC HII re- 
gion over the ones predicted from the isothermal, homo- 
geneous ionized gas model with two-density components, 
suggesting that dust radiation becomes significant in the 
continuum flux density. In order to evaluate the con- 
tribution of dust radiation in the UC HII regions, the 
terms involving the dust radiation B(T^) in Equations 
(1) and (2) are added to the B component model with 
the dust distribution as assumed. The flux density (Sd) 
of the dust radiation is determined using Equations (3) 
and (4) with an assumed dust temperature in the HII re- 
gion equal to the peak temperature in the Sgr B2 Main 
core, T d ~ 500 K, (|Lis fe Goldsmith II1990D and a power- 
law index for the dust opacity = 1.5 ([Goldsmith fe Lis I 
119901: iLis fc Goldsmith 1 1 1990( 1 . The observed brightness 
temperatures of few times 10 2 K for the continuum emis- 
sion from several central dense clumps also suggest a high 
dust temperature for the UC HII regions in Sgr B2 Main 
(also see the discussion in Section 5 for the high dust 
temperature). 

We show that, except for G, the dust radiation in the 
other seven regions makes a contribution to the observed 
continuum flux density at 0.86 mm comparable or larger 
than that from free- free emission (see Figured]). In gen- 
eral, the dust in each of the UC HII regions is optically 
thin, Td ~ 0.1, with a total gas mass < 11M Q . The dust 
optical depth and the gas mass (M to t) needed in each of 
the UC HII regions are summarized in Table 3. 

Finally, the models with two independent HII com- 



Zhao & Wright 



ponents mixed with dust involved a set of ten free pa- 
rameters, namely 2x(n e , T e , R, /y, and Tj). For F3, 
we have collected a total of nine measurements in line 
and continum flux densities from high-resolution observa- 
tions at centimeter-submillimeter wavelengths. Adding 
the two determined angular sizes for the components A 
and B, we have a total of eleven measurements to con- 
strain the model reasonably well. However, for the rest 
of the sources, the models need to be confined with fur- 
ther high-resolution observations of both continuum and 
line emissions at wavelengths from centimeters to sub- 
millimeters. 

4.5. Lyman continuum 

The hyper-compact HII components detected with 
H26a lines appear to require most of the flux of Lyman 
continuum photons from the ionizing stars to maintain 
their ionization. The H26a line flux from hyper-compact 
sources would be excellent tracers for newly formed mas- 
sive stars. Given an observed FWHM linewidth of 
AVfwhm that is dominated by Doppler broadening AVd 
in the H26a line profile, the ionizing photon flux required 
for the observed H26a line flux can be estimated 
from Equation (A8) that was generally derived for a ra- 
dio recombination line N + 1 — > N, 

(ifeH^)** *•-■<»> 

where as (T e ) is the total recombinati on coefficient to ex- 
cited levels (|Hummer fc Seaton II1963I ). and the function 
^26 provides a correction for the effects of both non-LTE 
and optical-depth. The values of \E , 26 are given in Table 
3 for each case. For the hyper-compact HII components 
(A) in Sgr B2 Main, * 2 6 is in the range ~ I (H26a-6) 
to - 2 (H26a-7). For the UC HII components (B) where 
the H26a line is optically thin, ~ b26, a correction 
for the factor due to the lower quantum-number levels 
that are underpopulated with respect to the LTE in the 
cases of lower electron density. The electrons at the level 
of N = 26 appear to be underpopulated by 20% to 30% 
with respect to the LTE for the UC HII regions in Sgr 
B2 Main. 

In addition, for H26a-1 (F3) assuming T c = 9,000 K, 
the value A^ym = 1-43 x 10 49 ph s _1 is inferred from the 
observed H26a line flux 280 Jykms -1 using Equation 
(5) with *26 = 1.3 while 2V Lym = 2.66 x 10 49 ph s" 1 is 
derived from the continuum flux densit y Sg = 2 - 9 Jy a t 
353.6 GHz using Equation (4) given by iWilcots I (| 1 994( 1 . 
The value determined from the observed continuum flux 
density surpasses the value determined from H26a line 
flux by 46%. Our detailed model fitting gives the free- 
free flux density of 1.61 Jy suggesting the fraction of dust 
contribution to the continuum is 44% at 353.6 GHz, in 
good agreement with above assessment. Considering op- 
tically thin of the line emission and nearly no attenuation 
and contamination by the dust at 353.6 GHz (0.85 mm), 
it appears to be an excellent way to determine the A^ym 
and the free-free flux density using the H26a line flux. 
Therefore, comparing the observed continuum flux den- 



sity with the free-free flux density determined from the 
H26a line flux, one can give a good assessment of the 
fraction of dust contribution to the continuum. 

4.6. Ionizing stars 

Based on the derived properties for both A and B 
components, we estimated the fluxes of the Lyman con- 
tinuum photons (A^Lym) from the newly formed ionizing 
stars. In each of the UC HII regions, the hyper-compact 
component (A) appears to require more flux in Lyman 
continuum photons than their larger but lower-density 
counterpart (component B) for the maintenance of the 
ionization (see Table 3). A total flux in Lyman contin- 
uum photons (A^f) is evaluated by the addition of the 
individual fluxes required for components A and B. As- 
suming that Aj^J 3 accounts for all the Lyman continuum 
photons from a single early-type, zero-age-main-sequence 
(ZAMS) star in each of the UC HII regions with no sig- 
nificant leakages, i.e. the UC HII regions are internally 
ionized (|De Pree et al. |[T998fl . the type of massive stars 
required for each of eight HII complexes with a hyper- 
compact H26a component is inferred on the basis of the 
stellar atmosphere model computed by iPanagial (|1973l ) 
and listed in Table 3. We note that the Lyman contin- 
uum photon fluxes from IPanagial (|1973H 's stellar atmo- 
sphere model were underestimated by 26% to 66% for 06 
to B0 type stars, respectively, as compared t o those com- 
puted from lVacca. Garmanv fc ShuU I dl996h 's model. In 
good agreement with lDe Pree et al. I (|1998f ) . we also find 
that at least an 06 star is required for the brightest H26a 
source (H26a-1) in F3 and an 08.5 star for H26a-2 in 
F4. 

The complex Fl, including H26a-3 (07), H26a-4 
(09.5) and H26a-6 (09), has been resolved by the VLA 
with resolution of .049"x 0.079" into at least seven com- 
pact components (|De Pree et al. I H9 98). The authors 
suggest that the complex requires a group of seven early- 
B and late-0 type stars (BO to 08.5) to maintain the ion- 
ization. The flux of Lyman continuum photons (A^ym) 
inferred from our analysis of the H26a line for the hyper- 
compact sources in F3 and Fl appear to be consider- 
ably greater than the values inferred from their 7 mm- 
counterparts based on the VLA continuum observations 
(|De Pree et al. IIT998D . The difference occurs because the 
hyper-compact components H26a-1 and H26a-3 in F3 
and Fl are optically thick in free-free emission at 7 mm, 
re ~ 2, inferred from the analysis above on the basis 
of H26a observations. The hyper-compact core appears 
to be deeply embedded in the optically thick region at 7 
mm. From the VLA flux densities at 7 mm, the flux of 
Lyman continuum photons has possibly been underesti- 
mated due to missing the contribution from the hyper- 
compact ionized core. Thus, an 07 star might be needed 
to maintain the ionization of hyper-compact H26a-3 in 
the HII complex of Fl. 

The above argument is also valid for the hyper- 
compact HII components H26a-5 (09), H26a-7 (09.5) 
and H26a-8 (08.5). These hyper-compact components 
appear to be embedded in the relatively isolated UC HII 
regions F10.37, F10.39 and G, respectively. In fact, the 
continuum optical depths inferred for the three hyper- 
compact HII components are greater than the values de- 
rived for the rest of hyper-compact HII components (see 
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Table 3). 

As discussed above (see Figure [3]), the H26a line from 
the three components appears to be enhanced by stim- 
ulated emission from continuum emission in the hyper- 
compact HII gas surrounding the newly formed O stars. 
The actual flux of the ionizing photons is reduced by a 
factor of 4 in H26a-7 as compared to an LTE source. 
Thus, the effect of stimulated line emission helps detec- 
tion of H26a line from a region ionized by an 09.5 among 
the eight H26a line sources in SgrB 2 Main. 

Given a hyper-compact HII region with T c = 1 x 10 4 K, 
AVd = 30 km s _1 without line stimulation, the 3cr (0.3 
Jy beam -1 ) detection limit for the H26a line imposed by 
the SMA data used in this paper gives a limit of N^ ym — 
7.1 x 10 47 ph s-\ or log(A Lym ) = 47.85, on the flux of 
ionizing photons, corresponding to an 09.5 ZAMS star. 
No significant detections of the H26a line with the SMA 
towards F2 are consistent with four B O stars inferred 
from the 7-mm continuum observations (|De Pree et al. I 
[19981) . 

4.7. Ages 

The ages of the newly formed O stars in Sgr B2 Main 
can be assessed by the dynamic time-scale of the larger 
HII components (B) as the time of sound wave traveling 
from the initial ionizatio n front at a radiu s close to the 
Str omgren r adius (R$) (IStromgrenl 119391) as s uggested 
by ISpitzer I (j!978l) and IGarav fe Lizanol (1 19991). U sing 
the Equation (3) given in Sh i. Zhao. fe Han I ([2010b) and 

the isothermal sound speed Cs = yx/ymn, we calcu- 
lated the dynamical age of the HII components (B) and 
listed the values in Table 3, ranging between the oldest 
one, 1.3xl0 3 y for F3 (H26a-1) and the youngest one, 
0.4xl0 3 y for F10.37 (H26a-7). If there is any exter- 
nal pressure produced by the swept-up ISM as the neb- 
ula expands, the sound crossing time will give an under- 
estimate of the age. The O type stars in Sgr B2 Main 
appear to be clustered around a thousand years ago. 

5. A CLUSTER OF PROTOSTELLAR CORES 

The dust emission from Sgr B2 Main has been resolved 
into at least twenty three components, showing that a 
cluster of protostcllar cores are present in this region in 
addition to the newly formed O and B type stars sug- 
gested by SMA observations of the H26a line and VLA 
observations of the radio continuum emission and H52a 
and H66a lines. From our interferometer observations, 
the brightness temperature (7b) of the continuum emis- 
sion from individual components can be determined un- 
der the Rayley- Jeans approximation, 



T b (i/) « 13.6A^ 



Sd 
Jy 



FWHM 



K. 



(6) 



where Sd is the dust flux density at wavelength, A from 
the region with the geometric mean of the angular size 
#fwhm in arcsec. Located near (0.13"±0.01"north) the 
brightest H26a line source in F3, Smm-1 has the high- 
est brightness temperature of 340 K. The contribution 
to the brightness temperatures in the continuum cores 
from the free-free continuum emission is less than 10%. 
Thus, the brightness temperature from dust emission is 
greater than 300 K in the center, suggesting that high 
dust temperature is present in Sgr B2 Main. 



From Equations(3) and (5), for a given dust tempera- 
ture Id, we can determine the optical depth (ra) for each 
of the dust cores from T^M: 



Td{v) = log 



Id 



Id - T h {v) 



(7) 



The dust temperature (Id) is affected by the increas- 
ing luminosity in the inner region (~0.1pc) of the Sgr 
B2 Main core where at least eight newly formed O stars 
are suggested by the H26a line sources. We determined 
the physical properties of the proto-stellar cores using a 
power-law distribution for the mean dust temperature Id 
as function of the radial distan ce r from the center. The 
distri bution of Td is mode l ed bvlScoville fc Kwan I (|1976[ ) 
and lWolfire fc Cassinelli 1 (fl987l) assuming that the total 
emissivity of the collection of grains at that tempera- 
ture equals the total radiative energy that is absorbed by 
the g rains (|Scoville fe Kwan|[T97llWolfire fc Cassinelli I 

T d = T in (—) (8) 



where T- m is the dust temperature at r; n = 0.01 pc, the 
radial distance between Smm-1 and H26a-1 (06 star) 
assuming a mean projection angle of 45°. Depending on 
the dust emissivity (Q^), the central stellar luminosity 
is insensitive to the exponent (3 of the dust emissivity 
power-law (r = 2/5 or 1/3 for /3 = 1 or 2, respectively). 

We used a linearly interpolated value of T = 0.37 
for j3 = 1.5 in the modeling. Assuming dust heated 
primarily by the central 06 star with a luminosity of 
2.5xl0 5 L o , T in (r in = 0.007pc) «400 K is inferred 
from the structur e of th e dust temperature derived by 
iScoville fc Kwan I (|1976| ). On the other hand, a larger 
value of T\ n (r\„ = 0.00 7 pc) ~ 630 K is calculated by 
iLis fc Goldsmith 1 (|1990l) in their detailed models assum- 
ing a total luminosity up to 2.3x10 7 Lq from the stars 
distributed in the core in Sgr B2 Main and exponent of 
/3 = 1.5 in the dust emissivity power-law. Because of 
high infrared opacity in the inner region of the core, the 
radiation from the outer region does not affect the dust 
temperature at the center. The observed brightness tem- 
perature, It, = 300 K from the dust in Smm-1 indicates 
the dust temperature Td ~ 470 K if Td w 1, which is con- 
sistent with the kineti c temperature of the absorbing gas 
towa r ds Sgr B2 Main ((Cernicharo. Go icoeche a. fc Pardo I 
I2006t lOin et al.ll2008[ ). We note that the dust sublima- 
tion radius i? S ub by heating from the central star can be 
determined by balancing the dust absorption of the UV 
and visible radiation from the star with the re-emission 
at IR from the dust grains at its sublimation tempera- 
ture T sub = 1800 K (IWolfire fc Cassinelli I fl986l) . For an 
06 star, i? su b = 2.9 x 10 15 cm, is about an order mag- 
nitude smaller than r m and in the zone between the two 
radial distances, the dust substantially cools down due to 
infrared radiation. In the following calculation, we take 
Tl n = 500 K in the power-law exponent of T = 0.37 for 
the dust temperature. Using Equation (8), we calculate 
Id in the range 500 K for Smm-1 near the 06 star to 170 
K for Smm-23 about 0.14 pc away. Using Equations (6), 
(7) and (8), we find that dust optical depth lies between 
0.3 and 3.0 at 0.86 mm for these protostellar cores, with 
unity for the brightest core Smm-1. Table 4 summarizes 
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the optical depths (r<i) 

From Equations (9), (10) and (11), we can calculate the 
total gas masses (Af), surface density (S) and molecular 
number density («h) in the individual protostellar cores 
assuming a mass ratio of gas-to-dust g = 100 and a mean 
mass per molecule, /i = 2.3, 



M = 12.0 



230GHz 



-0 



log 



T d - T b (i/) 



-Ckpc^FWHM 



A/,- 



! 



(9) 



£ = 94.3 



230GHz 



-/3 



log 



Ta 



T d - T h {v) 



g cm" 2 ,(10) 



and 



n H = 16.4 x 10 8 



v 



kpc"FWHM 



230GHz 
cm -3 . 



log 



T d - T h (v) 



(11) 



Table 4 lists the derived values for M, £ and nn in 
columns (6), (7) and (8). The model predicts a mass 143 
Mq for Smm-1, while the largest mass of 223 M Q is asso- 
ciated with Smm-4 near Fl (H26a-3). The enclosed gas 
masses (Af) of the individual dense cores range 5 to 223 
Mq. The molecular surface density is in the range 13 
to 150 gem -2 . The molecular number density is in the 
range Ixl0 8 -lxl0 9 cm -3 , which appears to be consis- 
te nt with the ambient den sity of 10 7 -10 8 cm -3 suggested 
by iDe Pree et aT~l (|1998| ) if ionized gas in the UC HII 
regions and the ambient m olecular gas are in pressure 
equilibrium (|Xie et al.lll996l ). 

In order to evaluate whether the individual dust cores 
are subject to inevitably undergoing gravitational col- 
lapse to form stars, we calculated the Bonnor-Ebert mass 
(A/be), the maximum mas s that a dense core remain s 
in hydrostatic equilibrium (jBonnor 195ft lEbert1ll957h . 
and free fall time scale (iff), 



Mre = 1.18CiG- 3/ V 1/2 



3tt 
2G~p 



(12) 



(13) 



where Cs, G and p are the sound speed, gravitational 
constant and mass density of the core. The derived values 
for A/be and tff are listed in columns 9 and 10 of Table 4. 
The masses of the protostellar cores in the central region 
of Sgr B2 Main are 1 to 2 orders of magnitude greater 
than their corresponding Bonnor-Ebert mass. Free-fall 
time scales of 1-3 x 10 3 y are inferred from our model 
calculation. 

6. DISCUSSION 

The kinematics of the Fl-4 HII complex and the dis- 
crete HII regions in Sgr B2 Main are shown by the radial 
velocity distribution determined from the first moment 
images of the H26a line in Figure |4] The observed kine- 
matics imply certain dynamical processes in the ionized 
gas associated with the newly born massive stars in the 
complex. Most of the H26a line emission arises from the 
high-density component (A) in an H26a source in Sgr B2 



Main, as shown with our model fitting (see Table 3). The 
high-density, hyper-compact components with a typical 
size of a few x 10 2 AU are candidates for rotating ionized 
disks. 

6.1. Candidate for rotating disk 

H26a-7 in F10.39, a weak stimulated source, shows a 
radial gradient, nearly north-to-south, with radial veloc- 
ity difference of 9 km s _1 or V T ~ 4.5 km s" 1 across the 
central 0.2" , or 2R — 1600 AU, which could indicate 
an edge-on rotating disk (see Figure [4} . If the velocity 
difference is due to Keplerian rotation of the ionized gas 
in a disk with an inclination angle i, then the enclosed 
mass is 



© 



V z 



km 



(lTO) Sin ~ 2i 







:18sirT 2 (i) M Q , 



(14) 



which is consistent with the mass of a late O-type star 
as the ionizing source. However, we note that the ob- 
served kinematics of the two sources can also be inter- 
preted as a rotating-expanding HII shell. The data are 
not adequate to discriminate between disk and rotating- 
expanding shell models. 

6.2. Confined HII ouflow and expanding ionized ring or 

disk 

F3 (H26a-1) shows a large velocity gradient northwest- 
to-southeast (see Figure 2]). The maximum radial veloc- 
ity difference is 50 km s -1 or U r « 25 km s -1 across 
a projected angular size of 0.7" or 2R = 5600 AU 
which is consistent with the angular size of the ring/disk 
seen in the VL A high-resolution (0.0 49" x 0.79") obser- 
vation at 7 mm (|De Pree et al. |[l998D . The ring/disk is 
nearly circular with a ratio of minor-to-major axis size 
^min/^maj > 0.5"/0.7", indicating an inclination angle 
i < 45°. If the velocity difference is due to Keplerian ro- 
tation of the ionized gas in the ring, then from Equation 
(14), an enclosed, dynamic mass 1.9 x 10 3 sin~ 2 (i) Af© 
can be inferred, which is at least an order magnitude 
greater than the inferred total mass of both the ZAMS 
type O star (~30Af Q ) in H26a-1 and the protostellar 
core Smm-1 (-34 Af ) in F3 (see Tables 3 and 4). On the 
other hand, the electron temperature of T e = 1.7 x 10 4 K, 
inferred for the large, lower-density component (B) from 
our model fitting, suggests an isothermal sound speed 
Cs = \J kT c /nift = 12 km s -1 which is about a factor 
of three smaller than the maximum expansion velocity 
Umax = U sin -1 (i) « 35 km s . However, the radiation 

pressure (P ra d = — -x) due to a luminosity L* from an 

3c7rit 2 

O type star would be comparable to the thermal pres- 
sure (Pth = 2n c kT c ) in the HII region with a radius of 
R, where the factor of 2 accounts for equal densities of 
electrons (n c ) and protons (rij), 



Prad =2.3 ^ 



L,. 



Pth 



\10 5 L C: 



.10 6 cm" 3 



R 



10 3 AU / 
V10 4 AV 



(15) 
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N F 



H26a-8 



o 

CD 

o 



6A 6.6 6 IB 1.0 ^ ^ 

o o 

CO 



F10.37 

H26a 



if 



a.- 

w 

o 
o 

H 
Q 



o 



LO 

d 




C/) 



1 0.5 -0.5 

RA offset (arcsec; J2000) 



F10.39 



1 0.8 0.6 0.4 0.2 



1.2 1 0.8 0.6 



FlG. 4. — The distribution of intensity weighted radial velocity with respect to a LSR velocity of 58 km s~ 1 is overlaid on the contour 
image of the continuum emission at 0.86 mm (sec Figure[T]l. The velocity image was constructed from the H26a line image cube convolved 
with a circular beam of 0.3" using the moment 1 algorithm in Miriad with 4<r intensity cutoff. The contours are (0.25,0.35,0.45,. . .,0.95) 
X 2.5 Jy beam -1 with a circular beam of 0.3". The cross signs mark the peak positions of the H26a line souces of 1, 2, 3, 4 and 6 in 
the central Fl-4 complex (the main figure) along with the sources of H26q-5, 7, and 8 in the discrete sources with large angular offsets 
from the center (the insets). Both the vertical and horizontal coordinates are the offsets (arcsec) from the registration center: R.A.(J2000) 
=17:47:20.135 and Decl.(J2000) =-28:23:04.53. 



where 



^ is the ratio of the radiation pressure to the 
thermal pressure. On the basis of the parameters de- 
rived for F3 (see Table 3), « 0.4 and 0.7 for the 

"th 

high-density component (A) and the lower- density com- 
ponent (B), respectively. On the other hand, if the radial 
velocity gradient is mainly due to expansion or outflow, 
the ram pressure of the expansion motions is puuV^ x . 
Then, the ratio of the ram pressure to the thermal pres- 
sure is, 



P 



= 6.1 x 10" 



th 



10 4 K 



\kms 



For the lower-density component (B) of F3. 



Pr: 



(16) 



0.5 if 



The derived ratios of ™ n and appear to 



Pth Pth 

be consistent. Both the thermal and radiation pressures 
play a critical role in the expansion of the ionized gas 
and may accelerate ionized outflow. However both the 
expansion and the plausible HII outflow are thought to be 



confi ned in the dense ambient medium (|De Pree et al. I 
fl99l) . If the UC HII regi on is in pressure equilibirium 
with the ambient medium ( Xie et al.lll996[ ). we have, 



2kn c T c t; = Tin (/iniHCTy + kT\ 



(17) 



where cry is the turbulent velocity, Tk is the kinetic tem- 
perature of the ambient gas, approximately equal to the 



dust temperature Td, and £ = 1 



P 



th 



considering the 



contribution from ram pressure. From Equation (17), we 
can find a minimum density nn required to confine the 
HII gas, 



n H >4 x 10 



? VlO 6 cm- 3 / V 



10 4 K 



AVi 



FWHM 



km 



where AVfwhm = \ (81n2) a 



cm 



(18) 



the FWHM 



line width of a molecular line from the surrounding 
medium. For Tk ~ Id — 500 K, the thermal width 
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is about 3 km s _1 . From the SMA measurements of 
the FWHM line widths for the sixteen kinematic absorp- 
tion features in H 2 CO(33o — 2 n2 ) and H 2 CO (3,^ — 2 2 o) 
towards the Fl-4 complex (Qi n et al.1 120081 ). we hnd a 
variance- weighted mean AVfwhm = 5±0.2 km s _1 , sug- 
gesting a dominant turbulent motion in the line broad- 
ening. Thus, from Equation (18) and the parame- 
ters derived for the HII regions (Table 3), a minimum 
7Jh ~ 2 x 10 7 cm -3 is needed to confine F3 (component 
B), which is a factor of a few to 10 less than the den- 
sities inferred for the protostellar cores but is in good 
agreement with the value of tih derived from the isother- 
mal model using the observed Tb = 14 K of the lowest 
contour in Figure [4] Our results agree with the pre- 
dicted range of the mo lecular nn between 10 7 — 10 s cm~ 3 
(|De Pree et al. Ifl99l . 

Fl (H26a-3, H26a-4, and H26a-6) shows compli- 
cated kinematics. The line flux from each of the three 
H26a sources suggests the presence of at least three O 
type stars in this regi on. From analysis of the high- 
resolution VLA image, iDe Pree et al. I (|1998| ) suggested 
a total of seven late O-type or early B-type stars in Fl. 
The observed kinematics suggest that complicated dy- 
namics of ioinzed gas are involved in the OB association. 

G (H26a-8) also shows a broad wing in its H26a line 
spectrum, which corresponds to a large radial velocity 
gradient (AV T ) northeast-to-southwest across the source 
with a size 0.6" as seen in the moment 1 image (see Fig- 
ure [4]) . Taking the parameters inferred from our model- 
ing (Table 3) and L*(08.5)« 5.4 x 10 4 L Q , with Equation 

(15), we find that * « 1 and « 0.6 for the hyper- 

compact component (A) and the larger component (B), 
respectively. The observed velocity gradient in the com- 
ponent B region appears to be mainly due to a bipolar 
ionized outflow accelerated in the region surrounding an 
08.5 star. The structure of kinematics for the hyper- 
compact component (A) with a weak stimulated source 
as suggested earlier was not resolved in these observa- 
tions. 

F10.37 (H26a-5) shows a broad H26a line profile with 
AVfwhm = 51 ± 6 km s _1 but no obvious radial veloc- 
ity gradient is seen in the moment 1 image (see Figure 
0]). For T e = 1 x 10 4 K, the thermal broadening gives 
AV^ WHM =21 km s _1 . Subtracting the thermal broad- 
ening, the residual line width AV F 1 W 1 HM =46 km s _1 re- 
quires large nonthermal motions. A bipolar ionized out- 
flow along the line-of-sight might be responsible for the 

large velocity dispersion ay = j= „ ^fwhm = 20 km 



s 1 in this region. 



V81n2 



6.3. Formation of massive protostellar cores 

In addition to the H26a sources, there are about two 
dozen protostellar cores with masses much greater than 
the corresponding Bonnor-Ebert masses, and a total 
molecular mass in the range 5 - 223 M Q . The observed 
filamentary structures in the continuum emission at the 
submillimeter wavelength suggest that the massive stars 
and protostellar cores were formed through fragamention 
during the initial collapse of the massive core in Sgr B2 



Main ()Qin et al.l 120081) . The inferred surface density of 
each protostellar core appears to be greater than unity; 
a minimum mass for forming a massive star in order 
to av oid further fragamentation (jKrumholz fc McKee I 
2008) . Through competitive accretion, each of the proto- 
stellar cores will grow by adding an accreted mass Mttf 
within the free-fall timescale t/ /. The accretion rate can 
be estimated for an isothermal flow in which the density, 
p(r) cx r~ 3 / 2 , and the enclosed mass within a radius r, 
M(r) cx r 3 / 2 , can be described by a power law (|Shu I 
I1977D . The accretion rate is linearly proportional to the 
accretion velocity Va C c and the enclosed mass Mr within 
r = R but inversely proportional to R, 



M = V A , 



dM(r) 



dr 

0.3A/QJ/- 1 



r=R 
V Ac 



kms 



Mr 

M Q 



-Y 

Auy 



(19) 



Based on our measurements 2R 
AU, M R - 1.2 x 10 3 M , and F Ac c 



V5" x 3" pa 3 x 10 4 



kT K 



1.3 



km s 3 for a kinetic temperature Tk 



/1TO H 

T d = 500 K 



we inferred M «0.03 M Q y . For a constant accretion 
rate of 0.03 M in a mean free-fall time-scale ?// ~ 
2 x 10 3 y about 60 M , or ~ 5% of the total mass 
increases in the protostellar cores. 

We have observed a snapshot of the process of form- 
ing a massive-star cluster. A few points regarding the 
formation of massive star cores can be drawn from our 
analysis of the SMA observations. First, Figure @] shows 
a noticeable angular offset of ~ 0.1" between the H26a 
sources (H26a-1, H26a-4 and H26a-6) from their contin- 
uum counterparts (Smm-1, Smm-4 and Smm-10). The 
relative angular offsets are significant (> 5<r), corre- 
sponding to a projected distance of 800 AU. If the H26a 
sources denote the locations of the new born O stars, 
the offsets of the protostellar cores (Smm-1, Smm-4 and 
Smm-10) from the ionizing stars might suggest that the 
local accretion centers are altered after the onset of the 
radiation from the O stars. A massive protostellar core 
in the immediate vicinity of a UC HII region has been 
also ob served in the W51e2 complex (|Shi. Zhao, fc Hani 
I2010albh . The authors use SMA high-resolution observa- 
tions to show that the protostellar core (W51e2-E) with 
a mass of 140 M Q , which is located 0.9" (4600 AU) from 
the UC HII region (W51e2-W), becomes the new accre- 
tion center and undergoes active star formation ejecting a 
powerful molecular outflow. Multiple pairs of O-type star 
and massive protostellar core found in Sgr B2 Main show 
a cluster version of the W51e2 system. Consequently, a 
considerable amount of mass must have existing in each 
of the paired protostellar cores before the O-type stars 
were born unless there were large accretion rates, ~0.1 
Mq s _1 at each of the massive protostellar cores in the 
past 1 x 10 3 years after the onset of radiation from the O 
stars. The protostellar cores and their companion mas- 
sive stars likely originate from fragmentation during the 
initial collapse of the cloud. After the primitive molecu- 
lar clumps formed, the massive ones undergo substantial 
gravitational collapse to form stars. Competitive accre- 
tion may occur in the process of forming a cluster of 
massive stars in a massive core like Sgr B2 Main. How- 
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ever, from our observations of this particular phase of 
massive star formation in a cluster, it is uncertain how 
much mass has been added to each of the protostellar 
cores through competitive accretion owing to the uncer- 
tainty in the time-scale of the process. On the basis of 
our analysis, accretion will add a small fraction of mass 
to each of the protostellar cores in a mean free-fall time- 
scale of 2x 10 3 y. 

7. CONCLUSION AND SUMMARY 

Using SMA archival data observed in 2007 for Sgr B2 
Main, we imaged the H26a line and continuum emis- 
sion at 0.86 mm with natural and uniform weighted syn- 
thesized beams of 0.40" x 0.28" and 0.36" x 0.22", re- 
spectively. Eight compact H26a line sources have been 
detected towards the central 5" (~ 0.2 pc). Using the 
non-LTE approach, we fitted the H26a line along with 
other available RRL and radio continuum data with 
two isothermal, homogeneous HII components: 1) high- 
density (a few xlO 7 cm -3 ), hyper-compact (several hun- 
dred AU) and 2) lower-density (a few xlO 6 cm -3 ), ultra- 
compact (a few thousand AU). The observed H26a line 
fluxes suggest that each of the eight H26a line emission 
regions must be powered by a Lyman continuum source 
from an O type star or a cluster of B-type stars. The 
typical age of the HII regions are a thousand years. The 
brightest H26a-1 source in F3 appears to be associated 



with an 06 star. The H26a line from H26a-7 (F10.39) 
appears to be enhanced by stimulated emission in the 
vicinity of an 09.5 star. 

About two dozen compact continuum emission cores 
are detected at 0.86 mm in the central 5" (0.2pc). Using 
a power-law distribution of dust temperature and ob- 
served brightness temperature, we determined the phys- 
ical properties for each of the protostellar cores with 
masses in the range of 5 to > 200 solar masses. The sur- 
face densities of the prostellar cores are in the range 13 
to 150 gem -2 , and molecular densities 1x10 s to lxlO 9 
cm -3 . We also calculated the Bonnor-Ebert mass (Mbb) 
for each core. The total molecular mass of each core 
appears to be much greater than the value of Mbe, sug- 
gesting that the cores undergo a substantial gravitational 
collapse to form massive stars with a typical free-fall time 
scale of a few thousand years. 

We are grateful to the anonymous referee for thought- 
ful comments and suggestions which were helpful for im- 
proving this paper. We thank Dr. Hui Shi for her assis- 
tance preparing Miriad scripts for the reduction of the 
SMA data. JHZ is grateful to the National Astronom- 
ical Observatories of China for hosting his visit during 
the course of writing this research paper. 



APPENDIX 

LYMAN CONTINUUM PHOTON FLUX 

In this Appendix, we show the detailed procedure to derive the Lyman continuum photon flux from the measurements 
of a hydrogen radio recombination line. 

Photoionizations and recombinations of hydrogen in HII regions 

For an HII region with electron density n c , the ionization equilibrium can be described by the balance between 
photoionizations and recombinations of electrons with the ions. The number of ionizing photons {N\) emitted by a 
star in unit time, or the Lyman continuum photon flux (A^ym), is equal to the number of recombinations (A^j, = 

^i? 3 n c 2 aB) to the excited levels within the Stromgren radius (R s ) in unit time ([Osterbrock 1119741 ): 

47T 

N LyirL = —R 3 s n 2 e a B (Al) 

where qb = C(-b(T c ) is the recombination coefficient in case B, which is a function of electron temperature (T e ) with a 
value a B = 2.58 x 10 -13 cm -3 s -1 at T c = 1 x 10 4 K (|Hummer fc Seaton I [19631) : and R 3 « w ^Cl ml D 2 Lfy. 

For an HII region located at distance D, with solid angle finiii path length L, and electron volume filling factor fy, 
Equation (Al) can be rewritten as, 

N Lym = n 2 e Lf v a B n ml D 2 (A2) 

Optical depth and flux density of hydrogen radio recombination lines 

For an isothermal- homogeneous HII region, the optical depths for hydrogen radio recombination lines (tl) and 
free-free continuum (tq), are proportional to the emission measure (EM = n 2 Lfy), and are given by [Shaver ( 1 9751 ): 

and 
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where r£ is the LTE line optical depth; &n and /3n are the population departure coefficients; n e is the electron density; 
AVd and AVp are the FWHM Doppler and pressure line widths in km s _1 , respectively. The correction factor a(y, T e ) 
is of order unity. 

On the other hand, if an HII region, with no dust in it, is optically thin and LTE (re << 1, t£ 1> an d 
&n = Asr = l)i the line flux density from Equation (1) in the main text becomes, 

2kis 2 

S*« si mi T e rZ (A5) 
&■ 

Normally, the Lyman continuum photon flux (A^L ym ) is derived fr om the emissio n measure determined from the 
free-free continuum flux density observed at a radio frequency, e.g. iWilcots 1 |l994). In the following, we show the 
derivation of the emission measure from a radio recombination line flux. Therefore, A^Lym can be determined directly 
from a radio recombination line flux. 

Emission measure from hydrogen radio recombination lines 

From Equations (A3) and (A5), neglecting the pressure broadening (AVd 3> AVp) the emission measure for an 
optically-thin, LTE gas can be expressed as, 

Jy kms~ 

For a general isothermal, homogeneous HII region, the emission measure can be determined by the following equation, 



^/ V = 5.66X10-WV 3^ K "».'. < A6 > 



"^=5.66x10-^-^0(3^) (^j)" 1 (|) 3/2 n„;,*;' (at) 

where *n is the ratio of the non-LTE line flux (Si,) to the LTE line flux in optically thin (S£[tZ « 1]) for the line 
transition N + 1 — > N, where N is the principal qauntum number, is a function of the line and continuum optical 
depths (tl and tc) and departure coefficients (&n and /?n), 



S£[t£ « 1] 



tl + t c 



1 - e -(^+rcA _ ( X _ 



L ) 



(A8) 



neglecting the term associated with the dust emission (S^) in Equation (1) of the main text. 



Lyman continuum photon flux from hydrogen radio recombination lines 

For an electron temperature (T e ) in an isothermal, homogeneous HII region, inserting n^Lfy of Equation (A7) into 
Equation (A2), we derived the Lyman continuum photon flux from the measurements with flux density (5l) and 
FWHM of Doppler brodening (AVd) for a hydrogen radio recombination line (N + 1 — > N) at frequency z^l, 

where \Pn corrects for the effects due to the optical depths and non-LTE distribution. Thus, from Equation (A8), 

e~ rc (l~e~ T L) 

~ 1 for optically thin, LTE gas; ~ &n(1 — ^tqPn) for optically thin, non-LTE gas; and *n ~ ; 

for the LTE gas becomes optically thick. 
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TABLE 1 

Log of observations and imaging 



Parameters 



H26a 



Continuum 



OBSERVATION: 
Date 

Array configuration 
Pointing R.A. (J2000) 
Pointing decl. (J2000) 

5u 

Bandwidth 
On-sourcc time 



2007 June 18 
Very extended (8 antennas) 
17:47:20.151 
-28:23:04.79 
348.664 GHz 
0.406 MHz 
2GHz(USB) & 2GHz(LSB) 
1.95 hr 
150-500 K 



CALIBRATORS: 
Flux density 
Phase 
Bandpass 



J1229+020, J1733-130. 



Callisto 
J1733-130 
J1743-038, 



J1751+096, J2015+371 



IMAGING: 
fo or A 
channel width 
FWHM beam 
r.m.s. 



353.623 GHz 
3 kms -1 
0.40"x0.28"(15°) 
0.1 Jy bm _1 chan - 1 



0.86 mm 

0.36"x0.22"(13°) 
8 mjy bin -1 
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TABLE 2 

SUBMILLIMETER PROPERTIES OF SCR B2 (MAIN) 



The H26ct Line Measurements 



Region 


Aa a 


A<5 a 


0Maj X 8 M in, PA 




Vlsr 


AVp W HM 


Sc 




(") 


(") 


(" X ", °) 


(Jy) 


(km s- 1 ) 


(km s- 1 ) 


(Jy) 


H26a-1 (F3) 


-0.17 


+0.03 


0.29x0.21, -25 


7.0 ±0.1 


64±1 


40± 1 


2.9±0.5 


H26a-2 (F4) 


+0.37 


+0.37 


U 


0.63+0.08 


66±2 


35±8 


0.41±0.10 


H26a-3 (F1A) 


-0.53 


+0.73 


U 


1.2 ±0.1 


48±2 


54±5 


1.8±0.4 


H26a-4 (FIB) 


-0.92 


+0.38 


U 


0.8 ±0.1 


73±1 


36±4 


1.0±0.4 


H26a-5 (F10.37) 


-0.00 


-1.18 


U 


0.62±0.08 


52±2 


51±6 


0.77±0.10 


H26o-6 (F10.303) 


-0.90 


+ 1.08 


U 


0.74 ±0.1 


57±2 


50±6 


0.91±0.24 


H26Q-7 (F10.39) 


+0.22 


-1.97 


u 


0.62±0.09 


69±2 


29±6 


0.35±0.09 


H26a-8 (G) 


+1.45 


+1.70 


u 


0.73 ±0.1 


53±2 


48±6 


0.2±0.1 



The Continuum Measurements at 0.8 mm 



Source 


Aa a 


A5 a 


SjVIaj X ^Min, PA 


Sp 


5*x 


A6> a 


T h 




(") 


(") 


(" X " °) 

0.80x0.40, +10 


(Jy beam x ) 


(Jy) 


(") 


(10 2 K) 


Smm-1 


-0.20 


+0.16 


2.6+0.2 


11±0.8 


0.26 


3.4 


Smm-2 


+0.78 


+0.09 


0.22 


0.4±0.1 


0.4±0.1 


0.79 


0.8 


Smm-3 


-0.39 


-0.72 


0.35x0.25, +3 


2.2±0.2 


2.2±0.2 


0.82 


2.5 


Smm-4 


-0.49 


+0.68 


0.83x0.30, +32 


2.0±0.3 


7.0±0.9 


0.84 


2.7 


Smm-5 


+1.00 


+0.41 


0.22 


0.4±0.1 


0.4±0.1 


1.08 


0.8 


Smm-6 


-0.92 


+0.59 


0.35x0.26, +21 


1.1±0.2 


2.2±0.3 


1.09 


2.4 


Smm-7 


+0.01 


-1.20 


0.61x0.31, -31 


0.9±0.1 


2.7±0.3 


1.20 


1.4 


Smm-8 


-0.68 


-1.09 


0.22 


0.8±0.2 


0.8±0.2 


1.28 


1.6 


Smm-9 


+0.58 


-1.17 


0.30x0.15, -75 


0.5±0.1 


0.6±0.1 


1.31 


2.0 


Smm-10 


-1.00 


+ 1.01 


0.38x0.22, -21 


1.3±0.2 


1.7±0.2 


1.42 


2.0 


Smm-11 


+1.03 


+ 1.04 


0.22 


0.3±0.1 


0.3±0.1 


1.46 


0.6 


Smm-1 2 


-0.08 


+1.59 


0.53x0.14, -82 


0.4±0.1 


0.8±0.2 


1.59 


1.1 


Smm-1 3 


-0.77 


-1.57 


0.53x0.45, -75 


1.2±0.2 


3.9±0.6 


1.75 


1.6 


Smm-14 


+0.33 


-1.74 


0.39x0.19, -48 


0.8±0.2 


1.3±0.3 


1.77 


1.7 


Smm-1 5 


-1.77 


-0.02 


0.45x0.09, +29 


0.4±0.2 


0.4±0.2 


1.77 


1.0 


Smm-16 


+0.83 


-1.67 


0.22 


0.4±0.1 


0.4±0.1 


1.86 


0.8 


Smm-1 7 


-0.89 


+ 1.64 


0.53x0.20, -13 


0.8±0.1 


1.6±0.2 


1.87 


1.5 


Smm-1 8 


-1.77 


+1.24 


0.22 


0.4±0.1 


0.4±0.1 


2.16 


0.8 


Smm-19 


+0.53 


-2.12 


0.22 


0.7±0.1 


0.7±0.1 


2.19 


1.4 


Smm-20 


-2.12 


+0.69 


0.22 


0.4±0.1 


0.4±0.1 


2.23 


0.8 


Smm-21 


+1.46 


+1.74 


0.22 


0.2±0.1 


0.2±0.1 


2.27 


0.4 


Smm-22 


-1.57 


+2.04 


0.30x0.19, +27 


0.5±0.2 


0.5±0.2 


2.57 


0.9 


Smm-23 


+0.56 


-2.57 


0.60x0.22, +57 


0.4±0.1 


1.0±0.1 


2.63 


0.7 



a Offsets from F3 at RA(J2000) = 17:47:20.178, Decl( J2000)=-28:23:04.67 



Table 3. Hyper- and Ultra- Compact HII cores in Sgr B2 (Main) 







H26a-1 


H26a-2 


H26a-3 


H26a-4 


H26a-5 


H26a-6 


H26a-7 


H26a-8 














(Fl-A) 


(Fl-B) 


(F10.37) 


(F10.303) 


(F10.39) 










A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


R I'm 3 ATH 




2.0 


4.0 


0.8 


3.2 


1.0 


4.0 


0.4 


2.7 


0.61 


2.4 


0.91 


2.8 


0.16 


1.2 


0.4 


2.8 


















HII GAS PROPERTIES 
















T. (K) 




9.0 


17 


10 


10 


9.5 


14 


5.0 


10 


5.0 


10 


5.0 


10 


7.0 


10 


6.5 


10 


n a (10 5 cm 3 ) 




42 


3.2 


50 


2.0 


52 


3.1 


85 


0.1 


52 


1.5 


30 


2.5 


330 


3.5 


110 


2.6 


IV 




0.1 


0.5 


0.1 


0.5 


0.1 


0.5 


0.1 


0.5 


0.1 


0.5 


0.1 


0.5 


0.1 


0.5 


0.1 


0.5 


b 2 6 




0.91 


0.82 


0.92 


0.73 


0.92 


0.79 


0.93 


0.67 


0.90 


0.72 


0.86 


0.74 


0.97 


0.76 


0.94 


0.74 


026 




-7.4 


-13.8 


-7.0 


-11.2 


-6.7 


-13.6 


-3.1 


-4.2 


-4.5 


-10.0 


-6.3 


-12.0 


-1.2 


-13.1 


-3.1 


-12.2 


T L (H26a) 




-0.60 


-0.008 


-0.30 


-0.012 


-0.33 


-0.015 


-1.1 


~0.0 


-0.61 


-0.003 


-0.36 


-0.011 


-1.3 


-0.02 


-0.71 


-0.014 


t c (0.86mm) 




0.024 


<0.001 


0.012 


<0.001 


0.019 


0.001 


0.040 


~0.0 


0.023 


<0.001 


0.011 


<0.001 


0.16 


<0.001 


0.050 


<0.001 


Ch26o 




94 


6 


73 


27 


71 


29 


99 


1 


94 


6 


81 


19 


95 


5 


81 


19 


*26 




1.3 


0.82 


1.1 


0.73 


1.1 


0.80 


1.7 


0.67 


1.3 


0.72 


1.1 


0.74 


2.0 


0.77 


1.4 


0.75 


r L (H66a) 




0.14 


-0.039 


0.058 


-0.085 


0.093 


-0.07 


0.71 


-0.001 


0.036 


-0.033 


0.019 


-0.068 


0.37 


-0.087 


0.32 


-0.073 


rc(13mm) 




9.5 


0.25 


4.8 


0.15 


7.3 


0.30 


4.2 


~0.0 


9.8 


0.065 


1.2 


0.052 


67 


0.044 


20 


0.23 


EM (10 10 cm" 6 


pc) 


1.6 


0.10 


0.97 


0.03 


1.4 


0.093 


1.4 


<0.001 


0.80 


0.013 


0.40 


0.042 


8.4 


0.035 


2.3 


0.045 


Mhii (10~ 2 M e 


) 


1.0 


3.3 


0.081 


1.0 


0.19 


3.2 


0.017 


0.030 


0.038 


0.33 


0.072 


0.85 


<0.005 


0.1 


0.022 


0.88 


U (pccm -2 ) 




74 


47 


34 


27 


45 


46 


24 


3 


27 


17 


28 


27 


24 


15 


29 


29 


log[jV L ym (ph S~ 


')] 


49.16 


48.31 


48.10 


47.81 


48.50 


48.34 


47.90 


44.98 


48.03 


47.19 


48.07 


47.82 


47.77 


47.02 


48.04 


47.85 






49.22 


48.28 


48.73 


47.90 


48.09 


48.26 


47.84 


48.26 


ZAMS-type 






06 




08.5 




07 


09.5 




09 


08.5 


09.5 


08.5 


*a g o (y) 






1.3 




1.1 




1.1 




0.9 




0.8 




0.9 


0.4 




0.9 
















DUST PROPERTIES (T d = 


500 K, 


= 1-5) 














T d (0.86mm) 






0.11 




0.03 


0.11 




0.17 




0.15 




0.12 


0.32 


<0.01 


M tot (M ) 






11 




1.9 




11 




7.5 




5.4 




5.9 


2.9 


<0.01 
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TABLE 4 

Properties of protostellar cores in Sgr B2 (Main) 













M 


2 






(10 J 




yp c ) 


/ir\3 ATT*\ 
AU ) 


(10 KJ 




( M <3) 


(g cm -2 ) 


cm ) 




Smm-l 


0.01 


4.5 


5.0 


1.1 


143 


55 


2.1 


2.9 


2.3 


Smm-2 


0.05 


1.8 


2.5 


0.4 


8 


20 


1.9 


1.0 


2.4 


Smm-3 


0.04 


2.4 


2.6 


2.6 


93 


130 


9.5 


0.5 


1.1 


Smm-4 


0.04 


4.4 


2.7 


1.9 


223 


92 


3.7 


0.9 


1.8 


Smm-5 


0.07 


1.8 


2.2 


0.4 


9 


22 


2.2 


0.9 


2.3 


Smm-6 


0.05 


2.4 


2.5 


3.0 


111 


150 


10.8 


0.4 


1.0 


Smm-7 


0.07 


3.5 


2.2 


1.0 


74 


48 


2.4 


0.8 


2.2 


Smm-8 


0.07 


1.8 


2.2 


1.3 


25 


63 


6.2 


0.5 


1.4 


Smm-9 


0.08 


1.7 


2.1 


0.9 


17 


46 


4.8 


0.5 


1.6 


Smm-10 


0.07 


2.3 


2.2 


2.3 


77 


112 


8.4 


0.4 


1.2 


Smm-11 


0.09 


1.8 


2.1 


0.3 


7 


17 


1.7 


0.9 


2.6 


Smm-12 


0.09 


2.2 


2.1 


0.7 


21 


35 


2.8 


0.7 


2.0 


Smm-13 


0.09 


3.9 


2.0 


1.6 


155 


79 


3.5 


0.6 


1.8 


Smm-14 


0.10 


2.2 


1.9 


2.1 


64 


105 


8.4 


0.4 


1.2 


Smm-15 


0.09 


1.6 


2.0 


0.6 


10 


32 


3.4 


0.6 


1.8 


Smm-16 


0.11 


1.8 


1.9 


0.6 


11 


27 


2.7 


0.6 


2.1 


Smm-17 


0.10 


2.6 


2.0 


1.4 


59 


68 


4.5 


0.5 


1.6 


Smm-18 


0.11 


1.8 


1.9 


0.6 


11 


28 


2.7 


0.6 


2.1 


Smm-19 


0.12 


1.8 


1.8 


1.5 


30 


75 


7.4 


0.3 


1.2 


Smm-20 


0.11 


1.8 


1.9 


0.6 


11 


28 


2.8 


0.6 


2.0 


Smm-21 


0.13 


1.8 


1.8 


0.3 


5 


13 


1.3 


0.8 


3.0 


Smm-22 


0.13 


1.9 


1.8 


0.7 


15 


33 


3.0 


0.5 


2.0 


Smm-23 


0.15 


2.9 


1.7 


0.6 


31 


28 


1.7 


0.6 


2.6 



